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ABSTRACT 

We show that accretion disks in several compact X-ray binaries with hydrogen- 
depleted donors are likely subject to a thermal ionization instability, unless they 
are strongly irradiated. These disks are particularly interesting in that their MHD- 
turbulent properties in the neutral phase may be quite different from those of standard, 
hydrogen-rich disks. 
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1. Introduction 

The vast majority of accretion disks around compact objects in close binary star systems 
are composed of hydrogen-rich material, as a result of Roche-lobe mass-transfer or wind-capture 
from a non-degenerate companion star (see Lewin, van Paradijs & van den Heuvel 1995 and 
Warner 1995 for reviews). It has long been known that such hydrogen-rich disks are subject 
to a thermal- viscous instability occurring when hydrogen becomes partially ionized, at central 
(mid-plane) temperatures T c ~ 10 K (or, equivalently, disk effective temperatures T e fj <^ 10 K; 
Meyer & Meyer-Hofmeister 1981; see, e.g., Cannizzo 1993a; Ludwig, Meyer-Hofmeister & Ritter 
1994; Hameury et al. 1998 for detailed calculations). 

There is also a well-established class of close binaries, the so-called AM CVn stars, in which 
the donor star is thought to be a helium white dwarf (e.g. Warner 1995). Helium disks in these 
systems are subject to the same thermal ionization instability as hydrogen-rich disks, as shown 
by Smak (1983) and subsequent workers (Cannizzo 1984; Tsugawa & Osaki 1997; El-Khoury & 
Wickramasinghe 2000). 
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Recently, the existence of a new class of close binaries, in which disks are not made of 
hydrogen or helium but instead metals such as and Ne (because of the nature of the white 
dwarf donors), has been proposed (Schulz et al. 2001; Juett, Psaltis & Chakrabarty 2001). In 
our study of the properties of supernova fallback disks (predominantly made of metals; Menou, 
Perna & Hernquist 2001a, b), we showed that disks composed of pure C or O are also subject to 
the thermal ionization instability. We argued that accretion disks of any metal content should be 
subject to this instability as well. 

In this Letter, we extend our work on metal disks to derive more precise global stability 
criteria. We then apply these results to the candidate metal disks in compact binaries, in an 
attempt to characterize their stability properties. 

2. Local Thermal Structure 

The thermal-viscous stability of gaseous, thin accretion disks around compact objects has 
been investigated extensively in the past. We follow a standard procedure here, similar to that 
employed to study the properties of fallback disks (Menou et al. 2001a). Because of the limited 
availability of metal opacities in the appropriate density and temperature ranges, our work is 
limited to C and O compositions only. The case of pure He is also considered for comparison with 
previous results. 

The Rosseland-mean opacities are taken from the OPAL database (Iglesias & Rogers 1996)Q. 
The opacities, at a specific mass density of 10 -6 g cm -3 , are shown as a function of temperature 
in Fig. |l], for the four compositions of interest: pure He (dash-dotted), pure C (dotted), pure O 
(dashed) and a C/O composition with a 40%-60% mass fraction (solid). The sudden opacity 
drop at temperatures T ^ 10 K, which corresponds to the recombination of the last free 
electron, is responsible in each case for the thermal ionization instability (Menou et al. 2001a). 
This temperature differs from composition to composition because of different ground state 
ionization potentials for each element. In particular, note that C dominates the opacity of the 
C/O composition at low temperatures (while O does at high temperatures). This is because C 
recombines at lower temperatures than O, and it implies that the stability properties of a disk 
with this C/O composition will be very similar to those of a disk made of pure C (see below). 

The disk thermal equilibria are found by equating the local viscous dissipation rate Q + to 
the radiative cooling rate Q~ , where Q~ is a function of the disk opacity. We neglect X-ray 
irradiation when calculating the thermal equilibria (see discussion below). The thermal equilibria 
are calculated from an extended grid of detailed models for the disk vertical structure, for different 
values of the surface density £ and the central temperature T c , at a given radius R from the central 
object of mass Mi (see Hameury et al. 1998 for details on the numerical technique employed). 



4 The OPAL database website is 
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Once a solution to the vertical structure (corresponding to a given S, T c and R) is known, the 
disk effective temperature T e $ is uniquely determined. The equation of state adopted for all the 
compositions is that of a perfect gas with a monoatomic adiabatic index 7 = 5/3. The lack of 
a detailed equation of state implies that vertical heat transport by convection, usually treated 
in the mixing length approximation, cannot be included in our calculations (see, e.g., Eq. [22] of 
Hameury et al. 1998). We note, however, that Gu, Vishniac & Cannizzo (2000) have argued for 
weak convection in thin accretion disks, because of the action of shearing and turbulent radial 
mixing. 

Figure § shows thermal equilibrium curves ("S-curves") thus obtained, in a surface density 
vs. effective temperature diagram. In each panel, the sections of the S-curves with negative slope 
are thermally and viscously unstable (see, e.g., Cannizzo 1993b and Lasota 2001 for reviews of 
the instability). For each composition, the thermal equilibrium solutions are shown for a different 
value of the viscosity parameter a, in order to guarantee that opacities are available in the range 
of density and temperature covered by the model. In the four panels of Fig. ^, the solid, dashed 
and dotted lines correspond to "S-curves" calculated at a radius R = 2 x 10 10 , 4 x 10 10 and 
6 x 10 10 cm, respectively, from a central object of mass M\ = 1.4M . (The calculation of only two 
S-curves for the He case was sufficient to test the consistency with previous results.) 

The instability occurs at a characteristic effective temperature T e ff iCr it (where the slope in the 
S-curve changes from positive to negative), which is a function of the disk composition and the 
radius of interest. When vertical transport of heat by convection is not included, as is the case here 
(because of the simple equation of state used), the value T e fj crit is also function of the viscosity 
parameter a (see below). As expected from the C-dominated opacities at low temperatures, the 
S-curves for the C/O composition are very similar to those of the pure C case. Our results for the 
He case are in good agreement with previous work on He disks (e.g. Smak 1983, Tsugawa & Osaki 
1997). 



3. Global Stability Properties 

The curves of thermal equilibrium can be used to derive a global stability criterion for the 
corresponding accretion disk. The relation M = 87rR 3 aT^ s /3GMi between the disk accretion 
rate and its effective temperature (e.g. Frank, King & Raine 1992) implies that a disk becomes 
thermally unstable if it locally accretes at a rate below the value corresponding to T e ff cr ; t (i2). 
The corresponding critical accretion rate is therefore obtained by fitting the turnover point of the 
"S-curve" at various radii for a given composition, as indicated by the dash-dotted line in Fig. [2]. 

For the four compositions of interest here, we find (after extrapolation to large a for the 
three metallic cases): 

He: M+ it (i?) ~ 5.9 x 10 16 mr ' 87 ^ 62 agf g s" 1 , (1) 
C: M+ it (R) ~ 1.3xl0 16 m r °- 7 Xo 23 «^ 4 g S ^, (2) 
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O : M+ it (R) ~ 5.0 x 10 16 mT/ " 68 ^ 05 «af g s _1 , (3) 
C/O: M+ it (R) ~ 1.2 x 10 16 mr°- 74 i*?6 21 «of g s" 1 , (4) 

where mi is the mass of the central object in solar units, R±q is the radius of interest in units of 
10 10 cm and cto.l = a/0.1. The scaling of Mi+ with mi was deduced by noting that mi enters the 
problem only via the square of the Keplerian angular frequency, Q 2 K oc mi/ ' R? , hence the direct 
proportionality to the scaling with radius (see Hameury et al. 1998 for the equations solved). We 
also have included in the above formulae the results of numerical explorations of the dependence 
of M^ it (R) on the viscosity parameter a (not explicitly shown in Fig. || for conciseness) 

These critical accretion rates do not differ strongly from the scaling M^- lt (R) ~ 
9.3 x 10 15 mj~°' 89 i?fg 68 g s" 1 found by Hameury et al. (1998) in the solar composition 
case (for a viscosity parameter a = 0.1). A disk of a given composition which satisfies everywhere 
M(R) > M^ it (R) is ionized and thermally (as well as viscously) stable. 

4. Applications 

In what follows, we apply the stability criteria derived in §3 to a few specific systems. In each 
case, we adopt a viscosity parameter a = 0.1 since there is empirical evidence from dwarf novae 
outbursts that this is the right order of magnitude for viscosity in a fully-ionized, MHD-turbulent 
thin disk (Smak 1999; Gammie & Menou 1998; Balbus & Hawley 1991; 1998). 

4.1. Metal-Rich Disk in 4U 1626-67 

Schulz et al. (2001) recently reported the Chandra detection of broad, double-peaked emission 
lines and strong photoelectric absorption edges of O and Ne in the low-mass X-ray binary pulsar 
4U 1626-67. Based on this, the authors argued that the disk in this ultra-compact system 
(with P orh = 42 min; Middleditch et al. 1981; Chakrabarty 1998) is probably fed by a C-O-Ne 
or O-Ne-Mg white dwarf of mass O.O2M . The short orbital period and small mass ratio in 
4U 1626-67 imply a disk outer radius ~ 2 x 10 10 cm (Schulz et al. 2001; Wang & Chakrabarty 
2001). 

Assuming a mass Mi = 1.4M for the neutron star and a viscosity parameter a = 0.1, we find 
that the mass-transfer rates below which the disk outermost regions in 4U 1626-67 become subject 



5 We have checked that, e.g. in the solar-composition case (for which a detailed equation of state was available), 
the dependence of M+ ;t on a nearly disappears when significant vertical heat transport by convection is allowed (see 
Fig. 2 of Menou et al. 2001a for a specific, quantitative comparison of stability in the presence and in the absence 
of convection). This result, which is confirmed by several previous vertical structure calculations (e.g. Cannizzo 
1993b; Hameury et al. 1998; Tsugawa & Osaki 1997), shows that our poor understanding of the role of convection 
for vertical heat transport in thin disks is a significant source of uncertainty for the derived stability criteria. 
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to the thermal ionization instability are ~ 4.7 x 10 16 g s" 1 and ~ 1.6 x 10 17 g s" 1 for pure C and 
pure O compositions, respectively. We note that the various compositions discussed by Schulz et 
al. (2001) for the white dwarf donor are dominated by either C or O, with some contribution from 
Ne. The above estimates of critical mass-transfer rates for stability should therefore be relevant 
in most cases because the disk stability properties tend to be dominated by the element with the 
lowest recombination temperature, as shown in §2 and 3 (Ne recombines first, followed by O and 
then C at lower temperatures; see e.g. Mihalas et al. 1990). However, if Mg makes a significant 
contribution to the disk composition (a possibility suggested by Schulz et al. for an O-Ne-Mg 
white dwarf donor), the critical mass transfer rate could be even lower because Mg recombines at 
lower temperatures than C (we were unable to address the Mg case in detail because opacities 
were not available). 

The above critical mass-transfer rates are between a factor a few and two orders of magnitude 
larger than the observationally- inferred value of ~ 10 16 g s _1 and the theoretically-inferred value 
of ~ 2 x 10 15 g s _1 discussed by Schultz et al. (2001) for 4U 1626-67. Our work therefore suggests 
that the accretion disk in this system should be subject to the thermal ionization instability^] 

4.2. Hydrogen-Poor Disks in Other Systems 

Juett et al. (2001) propose that four other low-mass X-ray binaries (4U 0614+091, 2S 0918- 
549, 4U 1543-624, 4U 1850-087) are also ultra-compact binaries with low-mass, hydrogen-depleted 
donors. Only 4U 1850-087 has a known orbital period of 20.6 min (clearly making it an 
ultra-compact system), but Juett et al. argue that the three other systems may also have orbital 
periods smaller than 1 hr. 

If this interpretation is correct (i.e. short orbital periods, low mass donors, helium- or 
metal-rich disks), we expect mass transfer driven by gravitational radiation at a rate similarily low 
as in 4U 1626-67 and disks equally subject to the thermal ionization instability.] 7 ] We note, however, 
that the knowledge of the binary orbital period is essential for the stability analysis because it 
largely determines the disk outer radius, which enters the stability criteria in Eqs. (l)-(4) with a 
large power. 

4.3. Helium-Rich Disk in 4U 1820-30 

The low-mass X-ray binary 4U 1820-30 is an ultra-compact system, with a 685 s orbital 
period, located in the globular cluster NGC 6624 (Stella, White & Priedhorsky 1987). The X-ray 

6 A similar conclusion is reached if the disk is assumed to have a solar composition. 

7 Again, the presence of a substantial amount of Mg in these disks could modify this conclusion. 
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burst properties of this system show that the donor star is a He white dwarf (Bildsten 1995). Given 
the short orbital period, the disk outer radius is ^ 10 10 cm and, according to Eq. (1), the critical 
mass-transfer rate below which the He disk would be subject to the thermal ionization instability 
is ^ 4.4 x 10 16 g s _1 . The inferred X-ray luminosity of the system is > 2 x 10 37 erg s _1 , however, 
implying an accretion rate in excess of 10 17 g s _1 (Bildsten 1995) and stability. This conclusion 
appears consistent with the inferred stability of the accretion rate in the system (Bildsten 1995), 
although the ~ 176 d low-amplitude luminosity cycle remains to be explained (Priedhorsky & 
Terrell 1984). 

5. Discussion 

We have calculated, as a function of radius, critical accretion rates below which disks become 
subject to the thermal ionization instability, for four specific hydrogen-poor compositions. We have 
applied these results to disks in several low-mass X-ray binaries with possibly Hydrogen-depleted 
donors and found that some of them could indeed be unstable. 

We caution, however, that our results are subject to a number of uncertainties. We have 
already mentioned the uncertainties associated with the magnitude of convection in thin disks 
(our results correspond to the inefficient convection limit). We have also shown that the exact 
composition of the disk material is important because the stability properties are in general 
dominated by the element with the lowest ground state ionization potential (e.g. Mg for an 
O-Ne-Mg white dwarf). Finally, we have neglected the stabilizing effects of disk irradiation in our 
calculations (van Paradijs 1996; Dubus et al. 1999). This is not unreasonable given that all the 
systems discussed are thought to be ultra-compact, but we note that irradiation may affect the 
stability of disks extending beyond 10 9 — 10 10 cm from the central neutron star, depending on 
their geometry (e.g. Dubus et al. 1999). 

The usual outcome of the thermal ionization instability in low-mass X-ray binaries and dwarf 
novae is the development of large amplitude outbursts and subsequent long periods of quiescence. 
However, this outcome crucially depends on the efficiency of angular momentum transport (or 
equivalently the value of the viscosity parameter a) in the disk when it becomes neutral (in the 
sense that a value of a significantly smaller in the neutral disk than in the ionized disk is required 
to produce large amplitude outbursts; Smak 1984).^] While there are good reasons for a to be 
significantly reduced in a neutral, hydrogen-rich disk because of MHD turbulence decay (Gammie 
& Menou 1998; Menou 2000; Fleming, Stone &: Hawley 2000), it may not be so in a hydrogen-poor 
disk. 



Several global numerical simulations have shown that, for similar values of a in the neutral and ionized phases, 
unstable disks experience only short, small amplitude luminosity variations (Smak 1984; Menou, Hameury & Stehle 
1999). 
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Menou (2002) emphasizes the importance of dissociative recombination by molecular hydrogen 
as a "sink" for free electrons and therefore a "source" of large resistivity in neutral, hydrogen-rich 
disks. In disks lacking hydrogen, the decay of MHD turbulence may be reduced or even absent, 
depending on the physical state and chemistry of the gas in the neutral phase. In our view, 
this is one of the most interesting aspects of the physics of hydrogen-poor disks, in the sense 
that a differential study of the disk stability and variability properties in hydrogen-poor and 
hydrogen-rich systems may reveal important clues on the nature of viscosity in quiescent, neutral 
disks. 
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log[T (K)] 

Fig. 1. — Rosseland-mean opacities as a function of temperature, at a mass density of 1CT 6 g cm -3 , 
for pure Helium (dash-dotted), pure Carbon (dotted), pure Oxygen (dashed) and a Carbon/Oxygen 
composition with a 40%-60% mass fraction (solid). The sudden opacity drop at T ^ 10 4 K 
corresponds to the recombination of the last available free electron. Note the similarity of the 
C/O opacity to the Carbon one at low temperatures and to the Oxygen one at high temperatures. 
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Fig. 2. — Examples of thermal equilibrium curves ("S-curves") for a thin accretion disk, shown in 
a surface density vs. effective temperature (S — T e g ) diagram. The curves correspond to a disk 
annulus located at 2 x 10 10 cm (solid line), 4 x 10 10 cm (dashed) or 6 x 10 10 cm (dotted) from 
a central neutron star of mass M\ = 1.4M . The four panels correspond to disks made of pure 
Helium (a), Carbon/Oxygen (40%-60% mass fraction; b), pure Carbon (c) and pure Oxygen (d). 
The value adopted for the viscosity parameter a and a fit to the critical effective temperature are 
shown in each panel. 



